In this paper an analytical solution is presented to solve Darcy-Brinkman equation for fully developed forced convection in a porous channel surrounded by two isoflux parallel plates, using homotopy perturbation method (HPM). The results show an excellent agreement between HPM and exact solution of the equation in various values of key parameters.

In this study, after discretization of the neutron diffusion equation and adjoint with high-order nodal expansion method in two dimensions and two energy groups, calculations with Momentum and Galerkin weighting functions for rectangular geometry (BIBLIS-2D) and hexagonal geometry (IAEA-2D) reactors are performed. The mean of relative power error for Momentum and Galerkin weighting functions was calculated in BIBLIS-2D reactor 0. 42% and 0. 62%, respectively, and for IAEA-2D reactor 4. 96% and 3. 52%, respectively. Regarding the results, it was concluded that in order to increase accuracy with the acceptable time of computing (4 Seconds for rectangular geometry and 28 seconds for hexagonal geometry with Intelr CoreTM i7-4510U Processor), the Momentum weighting function for rectangular geometry and the Galerkin weighting function for hexagonal geometry can be used to discretize equations without reducing the node size. Therefore, to increase the accuracy while maintaining the speed of calculations, without reducing the size of nodes, the appropriate weighting function can be used in discretization, which can be very useful in performing calculations of different transients.